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ABSTRACT 

According to our current cosmological model, galaxies like the Milky Way are expected 
to experience many mergers over their lifetimes. The most massive of the merging 
galaxies will be dragged towards the disc-plane, depositing stars and dark matter 
into an accreted disc structure. In this work, we utilize the chemo-dynamical template 
developed in Ruchti et al. to hunt for accreted stars. We apply the template to a sample 
of 4,675 stars in the third internal data release from the Gaia-ESO Spectroscopic 
Survey. We find a significant component of accreted halo stars, but find no evidence 
of an accreted disc component. This suggests that the Milky Way has had a rather 
quiescent merger history since its disc formed some 8-10 billion years ago and therefore 
possesses no significant dark matter disc. 

Key words: stars: abundances —- stars: kinematics and dynamics -— Galaxy: disc 
— Galaxy: formation — Galaxy: evolution — surveys 


* Based on data products from observations made with ESO 
Telescopes at the La Silla Paranal Observatory under pro¬ 


gramme ID 188.B-3002. These data products have been pro¬ 
cessed by the Cambridge Astronomy Survey Unit (CASU) at 


16 


© 0000 RAS 


2 Ruchti et al. 


1 INTRODUCTION 

In the standard ACDM cosmological model, galaxies build 
up their mass through the successive mergers of smaller 
galaxies. This process is known as hierarchical clustering 
(e.g. White & Rees 1978; Springel et al. 2006). It has long 
been a concern that the many mergers predicted by ACDM 
would destroy galactic discs (e.g. Stewart et al. 2008). In¬ 
deed, early simulations found excessive heating due to the 
bombardment of discs over cosmic time (e.g. Eke et al. 2000). 
More modern simulations have found, however, that the 
heating is dramatically reduced by moving to higher res¬ 
olution (Kazantzidis et al. 2008; Read et al. 2008) and if gas 
is included (Moster et al. 2010). The very latest cosmological 
simulations now resolve the most massive sites of star forma¬ 
tion ( ^ 100 pc) where most of the stellar feedback energy 
is generated; these are beginning to produce more realistic 
discs, for the first time (e.g. Guedes et al. 2011; Stinson et al. 
2013; Hopkins et al. 2014; Agertz & Kravtsov 2014). How¬ 
ever, even at this resolution potentially ad-hoc fine-tuning is 
required to suppress over-cooling in supernovae shocks (e.g. 
Stinson et al. 2006; Dalla Vecchia & Schaye 2008; Hopkins 
et al. 2014; Agertz & Kravtsov 2014; Kimm et al. 2015). 
We are closing-in on the issue of whether ACDM can be 
reconciled with the large abundance of disc galaxies in the 
Universe today, but the question remains open (e.g. Snyder 
et al. 2015; Crain et al. 2015). 

The Milky Way (MW) is a unique tool for understand¬ 
ing the formation and evolution of disc galaxies in the Uni¬ 
verse. We can study its stellar populations - which carry 
imprints of the history of the MW - in greater detail than 
in any other galaxy. By measuring the chemistry, ages, and 
dynamics of its individual stars we can trace their origins. 
This is known as Galactic Archaeology (Eggen et al. 1962; 
Ereeman & Bland-Hawthorn 2002). However, despite signif¬ 
icant progress due to the advent of large stellar surveys like 
the Sloan Digital Sky Survey (SDSS; Eisenstein et al. 2011), 
we still do not know the precise merger history of the Galaxy. 
This limits our ability to use the MW as a ‘rosetta stone’ for 
galaxy formation, since we do not know if its merger history 
is typical or rare. 

Minchev et al. (2014) argue that the relationship be¬ 
tween stellar velocity dispersion and o-abundance for stars 
from the RAdial Velocity Experiment (RAVE, Kordopatis 
et al. 2013) suggest early massive mergers, with a declin¬ 
ing merger intensity after disc formation some ^ 9 billion 
years ago; while Deason et al. (2013) argue that an observed 
‘break’ in the density profile of the MW’s stellar halo is also 
indicative of a quiet past, with few or no massive mergers 
over the past ~ 9 billion years. Both of these studies rely 
on model interpretations of the data but are sensitive to 
mergers over the full Galactic history. 

Another approach is to look for moving groups. Helmi 


the Institute of Astronomy, University of Cambridge, and by 
the FLAMES/UVES reduction team at INAF/Osservatorio As- 
trofisico di Arcetri. These data have been obtained from the Gaia- 
ESO Survey Data Archive, prepared and hosted by the Wide 
Field Astronomy Unit, Institute for Astronomy, University of Ed¬ 
inburgh, which is funded by the UK Science and Technology Fa¬ 
cilities Council, 
t email: greg@astro.lu.se 


et al. (2006) identified kinematic groups in the Galactic 
disc using the Geneva-Copenhagen Survey (Nordstrom et al. 
2004). At least some of these are due to resonant orbits in 
the disc (e.g. the Hercules stream; Bensby et al. 2007). How¬ 
ever, some could have extragalactic origins (Stonkute et al. 
2012, 2013; Zenoviene et al. 2014, 2015). The ages of the 
latter kinematic groups suggest they were accreted early on 
( > 8 Gyr ago). 

In this paper, we consider a new methodology that is 
most sensitive to massive ~ 1 : 10 mergers after disc forma¬ 
tion, but is less reliant on models: hunting for accreted stars 
in the MW disc (Ruchti et al. 2014, hereafter R14). Fur¬ 
thermore, since we look for the integral of all debris in the 
disc, we have greater sensitivity than other methods. The 
idea is already presented in detail in R14; here we briefly 
summarise the key points. 

Merging satellite galaxies will deposit stars, gas, and 
dark matter into the MW (e.g. Walker et al. 1996; Abadi 
et al. 2003; Read et al. 2008). For the most massive mergers, 
dynamical friction pulls the satellite towards the MW disc 
leading to an accreted ‘disc’ of stars and dark matter (e.g. 
Lake 1989; Read et al. 2008, 2009; Purcell et al. 2009; Ling 
et al. 2010; Pillepich et al. 2014). These stars - being born in 
shallower potential wells than the MW - should also show 
distinct chemistry (Helmi et al. 2006; Tolstoy et al. 2009; 
Kirby et al. 2011; Ruchti et al. 2014). 

If we can find these ‘accreted disc stars’, then they give 
us a powerful and direct probe of the merger history of the 
MW over the past ~ 9 billion years since the MW disc 
formed. The presence of such accreted disc stars encodes 
information about past massive mergers; its absence would 
suggest that such mergers did not occur. The accreted disc 
is also interesting because it contains dark matter - the so- 
called ‘dark disc’ (Lake 1989; Read et al. 2008, 2009). Such 
a dark disc, if present in the MW, has important implica¬ 
tions for the local dark matter density and velocity distribu¬ 
tion function. These in turn affect the expected dark matter 
signals in ‘direct’ and ‘indirect’ particle dark matter exper¬ 
iments (Bruch et al. 2009b; Bruch et al. 2009a; and for a 
review see Read 2014). 

In R14, we developed a new method to hunt for accreted 
disc stars. We applied our methodology to a small sample 
of about 286 stars with high-resolution abundances (Ruchti 
et al. 2011, 2013), finding little evidence for any accreted 
disc stars in the solar neighbourhood. We further demon¬ 
strated that this implies that the MW likely has a very light 
dark matter disc. However, the sample we used was small 
and nearby (< 1 kpc away), limiting the analysis to small- 
number statistics with strong observational and kinematical 
selection biases. We thus concluded that a larger, unbiased 
sample was needed to confirm our results. 

The Gaia-ESO Survey (Gilmore et al. 2012), which be¬ 
gan at the end of 2011, is an ongoing, five-year public spec¬ 
troscopic survey conducted with the FLAMES/GIRAFFE 
and UVES instruments at the Very Large Telescope. It will 
deliver high-quality spectra for ~ 100, 000 stars in the MW 
by 2017. The recently released internal third data release 
(iDR3) contains nearly 20, 000 spectra for stars in the MW 
field as well as globular and open clusters. Thus, the Gaia- 
ESO Survey presents an ideal data set with which to hunt 
for accreted stars in the MW. In this paper, we utilise this 
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extraordinary resource to expand our original analysis in 
R14. 

This paper is organised as follows. In Section 2, we 
briefly review the details of the chemo-dynamical template 
developed in R14 to identify accreted stars. In Section 3, we 
describe the spectroscopic data from the Gaia-ESO Survey, 
while the derivation of the physical and kinematic properties 
of the stars in our sample are detailed in Section 4. In Sec¬ 
tion 5, we search for substructure in the Gaia-ESO sample, 
identifying accreted stars. We then compute the mass and 
orbital parameters for the merging satellites in Section 6. 
Finally, in Sections 7 and 8, we discuss the implications of 
our results and give concluding remarks, respectively. 


2 BACKGROUND 

In this section, we briefly summarise the main aspects of our 
chemo-dynamical template. The reader is referred to R14 for 
more details. 


2.1 The kinematic template 

When satellite galaxies merge with the MW, they experience 
a dynamical friction force, given by: 

~ ( 1 ) 

v> 

where M^ax is the mass of the infalling satellite galaxy; v 
Is the deceleration due to dynamical friction; p is the back¬ 
ground density (i.e. stars, gas, dark matter etc.); C is some 
constant of proportionality; and v = |v| is the velocity of the 
satellite relative to the background (Chandrasekhar 1943; 
Read et al. 2006). From this equation, it can be determined 
that high-mass satellites will feel a larger force than low- 
mass satellites. Further, a larger force is expected for en¬ 
counters with the Galactic disc than the Galactic halo, since 
the former possesses a higher density of stars and gas. Thus, 
higher-mass satellites that encounter the disc will be “disc 
plane dragged” towards the Galactic plane over time and 
will deposit stars into an accreted disc-like structure. 

The accreted stars will be dynamically distinct from 
those stars that formed in the MW. This was shown in R14 
using simulations of a LMC-mass satellite galaxy under¬ 
going a collisionless merger with a MW-mass host galaxy. 
We considered varying orbital eccentricity and inclination 
as in Read et al. (2008). We investigated the location of 
the accreted and in situ stars in specific vertical angular 
momentum-energy space [JzjJc, Ez/Ec], where for a Galac¬ 
tic potential ^{R, Z) we have: 

Jz R-Ve Ez ^+^{R,Z) 

Jc Rc-Uc ’ Ez -h(R,0) ’ ^ 

and V& and Vz are the rotational and vertical velocity com¬ 
ponents with respect to the Galactic disc, respectively; 14 is 
the circular velocity, given by: 


Vc = 



( 3 ) 


and Rc is the radius of a planar circular orbit with the same 


specihc energy E as the star. This is found by minimising: 


^+-h(Rc,0)-F 


( 4 ) 


Note that a circular, in-plane orbit will have Jz/Jc = 
Ez/Ez = 1 . 

After the merger, the accreted population occupies re¬ 
gions of low- Jz/Jc and low-Fz/Rc, while the heated in situ 
stars primarily inhabit Jz/Jc > 0.8 and EzjEz > 0.95. How¬ 
ever, in R14 we found that the strongest kinematic differ¬ 
ential between heated and accreted stars is really seen pri¬ 
marily in Jz/Jc. For this reason, in this paper we will focus 
from here on solely on the Jz/Jc distribution. 

In Figure 1, we show Jz/Jc distributions from a simula¬ 
tion where a single LMG-like satellite merged on a prograde 
orbit inclined at 20° (left); 40° (middle) and 60° (right) to 
a Milky-Way like galaxy disc. The simulation was as out¬ 
lined in R14 and taken from the simulation suite already 
published in Read et al. (2008). In practice, we may expect 
several such mergers, each on a different orbit, leading to 
a broader or multi-peaked histogram rather than one clean 
near-Gaussian peak as in Figure 1. Here, we crudely estimate 
the effect of these multiple mergers by assuming identical 
mergers that linearly co-add stars. These are shown by the 
black (single); red (2-mergers) and blue (3-mergers) lines. 

The behaviour for LMG-mass mergers seen in Figure 1 
is distinct from lower mass mergers. As we move to lower 
mass satellites, we expect many more accretion events with 
a nearly Isotropic distribution (Read et al. 2008). Disc plane 
dragging by dynamical friction becomes less effective as M^at 
shrinks (see equation 1) and thus stars accreted from these 
low-mass satellite galaxies will have on average very low spe¬ 
cific angular momenta, peaking around Jz/Jc ~ 0.0. Such 
stars will contribute primarily to the canonical ‘stellar halo’, 
rather than the ‘accreted disc’. 

Based on these dynamical models, we concluded in 
R14 that ‘accreted disc stars’ should occupy regions where 
0.2 < Jz/Jc < 0.8, while stars accreted from small satellites 
typically end up as accreted halo stars (and/or streams) with 
Jz/Jc < 0.2. 

It is clear from Figure 1 that even a single LMC merger 
at 20° inclination shows up as a clear ‘bump’ in the Jz/Jc 
histogram. However, for even lower inclination angle, the 
merger debris will become increasingly buried in the tail 
of the Milky Way in-situ stars. As discussed in R14, the 
solution to this is to study also the distinct chemistry of 
accreted stars. 


2.2 The chemical template 

Stars born in small potential wells form from gas typically 
at lower surface density. This leads to a lower star formation 
efficiency and therefore lower alpha-to-iron ratios for a given 
metallicity bin (e.g. [Mg/Fe]^). This can be understood from 
simple chemical evolution models, but is also empirically ob¬ 
served in the surviving MW satellites (Matteucci & Brocato 
1990; Tolstoy et al. 2009, R14). It has also been shown that 
higher-mass satellites reach higher average metallicities (and 


^ The bracket notation is defined as [X/Y] = log(n(X)/n(Y))* — 
log(n(X)/n(Y ))0 for any elements X and Y. 
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Figure 1. ‘Gaia-ESO volume’ (6 < i? < lOkpc; \z\ < 2kpc) heated and accreted stars taken from a collisionless simulation of an LMC- 
like prograde merger with a MW-mass disc galaxy (Read et al. 2008). The initial orbital inclination of the merging satellite was chosen 
to be 20° (left); 40° (middle) and 60° (right). The panels show histograms of JzjJc (the ratio of the star specific angular momentum 
about the disc axis to that expected for a circular orbit at the guiding radius) for all stars assuming a single LMC-like merger (black) two 
identical such mergers where the stars are assumed to co-add linearly (red); and 3 such mergers (blue). The distributions are normalized 
such that their total area is equal to unity; they are all well-fit by a double Gaussian. In all cases, the accreted stars show up as a clear 
bump at lower Jzj Jc. than the in-situ star distribution. For very low inclination and/or low eccentricity mergers, however, the kinematics 
of the accreted and in-situ stars can overlap. For this reason, chemistry becomes important since the metallicity and abundance of stars 
born in a satellite should differ from those born in-situ in the Milky Way disc (see §2.2). 


possibly higher alpha-to-iron ratios for the largest masses) 
than lower-mass satellites (e.g. Lanfranchi et al. 2006; Kirby 
et al. 2013). In R14, it was shown that stars in present- 
day surviving dwarf galaxies, as well as the LMC, typically 
have [Mg/Fe] ^ 0.3 above [Fe/H] ^ — 1.5 (see their Fig. 1). 
We therefore expect that approximately LMC-mass mergers 
that will be significantly disc plane dragged, will also have 
[Mg/Fe] ^ 0.3 at metalliclties above [Fe/H] ~ —1.5. 


The combination of chemistry and kinematics is re¬ 
quired to unambiguously detect accreted stars in the disc. 
Stars can also be born in-situ at low efficiency and low metal¬ 
licity in the outer disc, migrating later on into the solar 
neighbourhood. Similarly, accreted stars from satellites on 
particularly planar, circular, orbits could masquerade kine¬ 
matically as in-situ stars (R14). For these reasons, it is the 
combination of ‘hot’ stellar kinematics and low alpha abun¬ 
dance that provides the smoking-gun evidence for accreted 
stars that we searched for in R14 and that we will search for 
in this paper. 


In R14, we applied our analysis to a small sample of 
286 stars, and restricted our search for accreted stars to only 
regions where [Mg/Fe] < 0.3. However, with a large enough 
sample, we can search for sub-structure in the distribution 
of Jz/Jc over all the metallicity-abundance space. In the 
present work, we use a sample from the Gaia-ESO Survey 
that is more than ten times the size of the sample in R14 to 
search for accreted stellar components in the MW. 


3 OBSERVATIONAL DATA 

In this paper, we make use of the iDR3 data from the Gaia- 
ESO Survey for 7, 723 MW field stars with the GIRAFFE 
spectrograph, which have been observed with both the HRIO 
(i? ~ 19,800) and HR21 (R ~ 16,200) settings. All targets 
were selected according to their colour and magnitude, using 
the Visible and Infrared Survey Telescope for Astronomy 
(VISTA) photometry. The selection function of the Gaia- 
ESO Survey will be presented in a future paper (Gilmore 
et ah, in prep.), while some details can be found in Recio- 
Blanco et al. (2014). 

Radial velocities are measured from the GIRAFFE 
spectra, with typical uncertainties of ~ 0.3 km s“^ (see Ko- 
posov et al. in prep., for more details). The stellar effective 
temperature (Teff), surface gravity (log g), and metallic¬ 
ity ([Fe/H]) were derived by the Gaia-ESO Survey work¬ 
ing group in charge of the GIRAFFE spectrum analysis for 
FGK stars (WGIO). The details of this analysis will be pre¬ 
sented in Recio-Blanco et al. (in prep.). Briefly, the spectra 
were analysed by three research groups using different spec¬ 
tral analysis codes: MATISSE (Recio-Blanco et al. 2006), 
FERRE (Allende Prieto et al. 2006), and SME (Valenti & 
Piskunov 1996). 

All analyses were performed using the ID LTE 
spherically-symmetric (log g < 3.5) and plane-parallel 
(log g ^ 3.5) MARCS model atmospheres (Gustafsson et al. 
2008) and a common atomic (and molecular) line list (Heiter 
et al. in prep.). The results from the three different groups 
were corrected for method-to-method offsets using multi¬ 
linear transformations as a function of Teff and log g, and 
then averaged to give a final estimate of the stellar param¬ 
eters for each spectrum. The stellar parameters were also 
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Figure 2. Location of our sample, consisting of 4, 675 stars, in the 
Hertzsprung-Russell diagram. The stars are color-coded according 
to their measured [Fe/H]. 


calibrated with respect to several benchmark stars (Jofre 
et al. 2014). Subsequently, individual abundances were es¬ 
timated from the spectra and averaged amongst the three 
methods. 

For our analysis, we make particular use of [Fe/H] and 
the [Mg/Fe] ratio. We thus utilized the measured Fe and 
Mg abundances from Fel and Mgl lines in the GIRAFFE 
spectra, respectively. Moreover, we averaged the results for 
several measurements of the solar spectrum to compute the 
solar abundances. We found average values of A(Fe )0 = 7.51 
and A(Mg )0 = 7.60. These values are very similar to the 
solar abundances measured by Asplund et al. (2009), which 
are the solar abundances used in R14. We therefore used 
these solar values when computing [Fe/H] and [Mg/Fe] in 
order to place the data on similar scales to R14. 

Since the quality of the final stellar parameters de¬ 
pends on the signal-to-noise ratio (SNR) of the spectrum, 
we selected only stars with SNR 15 in both HRIO and 
HR21. Additionally, to reduce possible abundance scatter¬ 
ing in our analysis, we further restricted our sample to only 
those stars with [Mg/Fe] uncertainties of ^[Mg/Fe] < 0.2. 
We further removed four stars that did not have a mea¬ 
sured Mg abundance. Finally, we removed any stars that 
do not have 2MASS or VISTA magnitude estimates and 
UCAC4 (Zacharias et al. 2013) or PPMXL (Roeser et al. 
2010 ) proper motions, since these are important for deriv¬ 
ing the kinematics of the stars (see §4.3). 

Our final sample consists of 4, 675 stars, with average 
uncertainties in their final stellar parameter estimates of 
(St,,,) ~ 68 ± 75 K, (Jiog g) - 0.13 ± 0.11, and (<5[Fe/H]> ~ 
0.06 ± 0.03, while (5[Mg/Fe] ~ 0.05 ± 0.02. These uncertain¬ 
ties were estimated from the method-to-method scatter in 
the combination process. Figure 2 shows the Hertzsprung- 
Russell (H-R) diagram for our hnal data sample. 


4 THE DYNAMICAL ANALYSIS 

4.1 VISTA vs. 2MASS magnitudes 

For most stars, both VISTA and 2MASS magnitudes were 
available in the Gaia-ESO internal data archive. The VISTA 
photometry typically has a higher precision than 2MASS, es¬ 
pecially for dim stars (mean uncertainty in Ks is ~ 0.03 mag 
versus ~ 0.06 mag, respectively). Thus, we adopted the 
VISTA Jv and Ks,v magnitudes when available for a given 
star. There were four cases in which no VISTA magnitude 
was listed. For these, we adopted the 2MASS magnitude. 

Our derivation of the distance to the stars requires a 
2MASS 7 ^ 3,2 magnitude (see §4.2). Thus, for stars in which 
we adopt the VISTA photometry, we converted the VISTA 
Jv and Ks,v magnitudes to the 2MASS J 2 and As, 2 magni¬ 
tudes using: 

As ,2 = As,v — 0.01081 (Jv — As,v) (5) 

J 2 = Jv - 0.07027 (Jv - As,v). (6) 

These equations were derived from those that were deter¬ 
mined to account for zero-point offsets between the VISTA 
and 2MASS magnitudes^. 

Although we adopt a combination of VISTA and 
2MASS photometry in order to maximise the size of our 
sample, our overall results do not change if we adopt one or 
the other, exclusively. 


4.2 Absolute magnitudes and ages 

Both the age and absolute magnitudes Mks 2 of 

each star in our sample were inferred from the estimated 
stellar parameters (Teft, log g, and [Fe/H]) and a library 
of evolutionary tracks by using BeSPP (Bellaterra Stellar 
Parameter Pipeline). A detailed description of the Bayesian 
method implemented in BeSPP can be found in Serenelli 
et al. (2013). 

In this work, we use a newly computed grid of stellar 
models using Garstec (Weiss & Schlattl 2008). While most 
of the physical input in the models is the same as described 
in Serenelli et al. (2013) there is one important exception. 
For the new grid of models, certain amount of extra-mixing 
below the convective envelopes has been added following a 
phenomenological recipe similar to that given by Vanden- 
Berg et al. (2012). Further details will be discussed elsewhere 
(Serenelli 2015, in prep.). 

A solar model has been calibrated to fix the mixing 
length parameter for all the grid models (omut = 1.798). 
The initial composition of the solar model is Zq = 0.01940 
and Vq = 0.2712. For the whole grid we choose the GN93 
solar composition (Grevesse & Noels 1993) because it leads 
to solar interior models that are in agreement with helio- 
seismic data. As shown by Ghristensen-Dalsgaard (2009), 
this choice leads to a correct estimation of the solar age, 
in contrast to models using the more recent AGSS09 solar 
composition (Asplund et al. 2009) which overestimate the 
solar age by about 10 %. 

The grid of stellar models spans the mass range 0.6 ^ 

^ For more details, see: http://casu.ast.cam.ac.uk/surveys- 
projects / vista / technical / photometric-properties 


© 0000 RAS, MNRAS 000, 000-000 



6 Ruchti et al. 


M/Mq ^ 3.0 with steps of AM = 0.01 Mq and the metallic- 
ity range —3.0 ^ [Fe/H] ^ 0.6 with steps of A[Fe/H] = 0.1 
for [Fe/H] < 0 and A[Fe/H] = 0.05 for [Fe/H] > 0. The 
initial composition of models is fixed by assuming a con¬ 
stant Ay/AZ = 1.17 as determined by anchoring this rela¬ 
tion to the solar model and the helium abundance Tsebn = 
0.2485 determined from standard Big Bang Nucleosynthesis 
(Steigman 2007). 

As it has been originally discussed by Salaris et al. 
(1993), a-enhancement can be mimicked in solar-scaled stel¬ 
lar tracks by modifying [Fe/H] by an amount proportional 
to the a-enhancement. Based on a subset of stellar tracks 
including a-enhancement, we find the following relation for 
the adjusted [Fe/H] 

[Fe/H]adjusted = [Fe/H] -b 0.22 [Mg/Fe]/0.4, (7) 

where ]Mg/Fe] is used as a proxy for the Q-enhancement 
of the star. Therefore, in BeSPP we effectively account for 
Q-enhancement by fitting T^s, log g, and ]Fe/H]adjusted. 

In Serenelli et al. (2013), we showed that the stellar 
mass and age cannot be determined for stars with log g < 
3.5. This is due to the degeneracy of the stellar evolutionary 
tracks in the Teff — log g plane at these surface gravities. 
Therefore, in this work we only present ages for stars with 
log g > 3.5. 

From these fits, we obtain a full probability distribution 
function (PDF) for the absolute Mk^ 2 magnitude and age 
of a given star. These PDFs are not Gaussian for all stars. 
We can thus use these PDFs to determine snbsequent PDFs 
for the kinematic and orbital information for the stars. 


4.3 Distances, kinematics, and orbital parameters 

The distances, kinematics, and orbital parameters of the 
stars were derived using a Monte Carlo approach. First, 
we constructed PDFs for the J 2 and Ks ,2 magnitudes (de¬ 
scribed in §4.1), proper motions, and radial velocities (from 
the Gaia-ESO Survey) assuming Gaussian distributions with 
sigma equal to the uncertainty in the given quantity. For a 
single star in our sample, we next created 3000 Monte Carlo 
sets of these parameters drawn from these PDFs, as well as 
the (possibly non-Gaussian) PDFs for the absolute Mxg 2 
magnitude. We then derived the distance, the 3D space mo¬ 
tions, and orbital parameters for each Monte Carlo param¬ 
eter set. Each of these derivations thus has an associated 
PDF. Below we give further details for each of the parame¬ 
ters. 


4 . 3.1 Distance 

The distance was determined following the iterative proce¬ 
dure described in Ruchti et al. (2011). This method uses 
Ks ,2 and Mk^ 2 derived an initial distance estimate 
using the distance modulus. The reddening is then esti¬ 
mated from the Schlegel et al. (1998) dust maps. Values 
oi E{B — V) > 0.1 were corrected according to the for¬ 
mula given in Bonifacio et al. (2000): E{B — V)corrected = 
0.035+0.65*77(5 —V)schiegei. Finally, we assume that for dis¬ 
tances within the dust layer (assumed to have a scale-height 
of h — 125 pc), the reddening at a distance D and Galactic 
longitude b will be reduced by a factor 1 —exp(—]Dsinfe]//i). 


This reduction only affects those stars that are relatively 
nearby and lie close to the Galactic plane. The E{B — V) 
reddening was then converted to the A's _2 extinction coef¬ 
ficient, given by Ak^ 2 ~ 0.366 * E{B — V) (Casagrande & 
VandenBerg 2014). The mean value of Aks 2 1®®® than 

0.1 mag, with an average uncertainty of 0.1 mag. This 
corresponds to a .5], 5% uncertainty in the distance. The 
distance estimate and reddening were iterated until there 
was a difference in distance of 1% between two consecutive 
iterations. 


4 . 3.2 Kinematics and orbital parameters 

We next combined the distance with the proper motions and 
radial velocities to derive the full 3D space motions in Galac- 
tocentric cylindrical coordinates {Vn, Ve, Vz). We adopted 
UGAC4 proper motions for the stars. However, when un¬ 
available (this is especially true for the more distant stars), 
we used the PPMXL proper motions. Although there are 
possible systematics between the two sets of proper motions 
(see Zacharias et al. 2013), the overall results of our analysis 
do not change if we adopt one or the other. 

The average uncertainties in the proper motions are 
~ 5 mas yr“^ and ~ 7 mas yr“^ for UCAG4 and PPMXL, 
respectively. The velocities were corrected to the solar ve¬ 
locity with respect to the local standard of rest (LSR), 
given by (Uq,Vq,Wq) = (14.00,12.24,7.25) km s“^ from 
Schonrich et al. (2010) and Schonrich (2012). An additional 
Vlsr = 220 km s“^ was added to put the space motions in 
the Galactic rest frame. 

We next computed the orbital parameters of the stars 
following the methodology described in R14, which assumes 
a three component Galactic potential: a Hernquist (1990) 
bulge, a Miyamoto & Nagai (1975) disc, and a Navarro et al. 
(1996) dark matter halo. For these calculations, we assumed 
the Sun sits at Rq = 8 kpc. We computed the specific an¬ 
gular momentum (Jz/Jc), as defined in §2. 

The advantage of using Jz/Jc instead of other orbital 
parameters, such as the orbital eccentricity, is that it is less 
dependent on the chosen Galactic potential. For example, 
we have adopted the masses and constants for each poten¬ 
tial from Gomez et al. (2010). However, it is possible that 
the bulge mass was over-estimated (e.g. Laurikainen et al. 
2010). If we lower the mass of the bulge by half, this leads 
to differences of less than 0.05 in Jz/Jc. 

4 . 3.3 Final PDFs and uncertainty estimates 

Following the procedures described in the previons sections, 
we then combined the estimates for all 3000 Monte Garlo pa¬ 
rameters sets to construct the PDFs for the distance, space 
motions, and Jz/Jc for each star in our sample. The best 
estimate for each quantity was taken to be that at which 
the PDF is maximised. The lower and upper uncertainties 
were adopted to be the values a the boundary of the 68.3% 
confidence interval snrrounding the value at the maximum. 

This resulted in typical upper and low age uncertainties 
of ~ 4 Gyr and ~ 1 Gyr, respectively. The typical uncer¬ 
tainty in the distance was found to be 10%. This, com¬ 
bined with uncertainties in the proper motions and radial 
velocities, led to an the uncertainty in the space velocities 
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that is typically < 50km s~^. All of this together leads to 
mean upper and lower uncertainties in Jz/Jc of ~ 0.11 and 
~ 0.08, respectively. 

Larger uncertainties in the distance and proper motions 
often lead to increased uncertainty in Jz/Jc. For many stars, 
however, the PDFs for the distance are not necessarily Gaus¬ 
sian. This is also the case for the age PDFs. Further, corre¬ 
lations among the different uncertainties are very non-linear 
and complex. For example, the distance and proper motions 
are determined independently. Thus, larger uncertainties in 
the distance do not necessarily correspond to larger uncer¬ 
tainties in the proper motions. It is therefore important to 
retain the information in the derived PDFs for Jz/Jc and 
ages during our analysis. 


4.4 Aside: the Gaia-ESO selection function 

In R14, our sample was quite small and kinematically bi¬ 
ased due to the selection function in the Ruchti et al. (2011) 
study. Our current sample, however, is more than ten times 
the size of the R14 sample. Moreover, there is no kinematic 
selection for the MW fields in the Gaia-ESO Survey. Instead, 
the stars are selected according to their colour and magni¬ 
tude. This can cause some bias in the age—metallicity rela¬ 
tion of the stars, in the sense that young, metal-poor stars 
and old, metal-rich stars are under-sampled (e.g. Bergemann 
et al. 2014). Further, the selection means that the stars in 
our sample probe different volumes, with dwarfs and giants 
probing smaller and larger volumes, respectively. This means 
that our sample is not volume complete, and so our sample 
may over- (or under-) populate stars from different regions 
of the MW. This can cause a bias in the overall kinematic 
distributions and number density of stars in relation to the 
distance from the Sun. This is due to the fact that the typ¬ 
ical motions of the stars in the MW differ, depending on 
their position in the Galaxy. 

Given our sample is not volume-complete, we cannot 
easily quantify the relative number density of accreted stars 
versus m situ stars within a specific volume. However, as 
shown in R14, stars identified as accreted stars in our chemo- 
dynamical analysis are all genuine accreted stars, even for 
a more complicated selection function, which includes kine¬ 
matic biases. Since there is no kinematic selection in the 
Gaia-ESO Survey, it will have a small effect on our present 
analysis. We therefore can study the dynamical distributions 
to identify accreted stars in our sample without any under¬ 
lying corrections to the models presented in §2. 


5 HUNTING FOR SUBSTRUCTURE 

5.1 Abundance trends of the low Jz/Jc stars 

The size of our present data set from the Gaia-ESO Survey 
allows us to search for dynamical substructure throughout 
the [Mg/Fe] — [Fe/H] space. As an initial exploration, we 
first looked at the abundance trends of the stars that have 
orbits consistent with possible accretion, i.e. Jz/Jc < 0.5 
(see Figure 1). 

In Figure 3, we plot our sample in the [Mg/Fe] — [Fe/H] 
plane, highlighting those stars with Jz/Jc < 0 (blue) and 
0 ^ Jz/Jc < 0.5 (red). We further separate the sample into 
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Figure 3. Location of stars in the [Mg/Fe] — [Fe/H] plane, high¬ 
lighting stars with Jz/Jc < 0 (blue triangles) and those with 
0 Jz/Jc < 0.5 (red circles), (a) All stars in our sample, (b) Only 
stars with log g /S 3.5, defined as sub-giants/dwarfs, (c) Only 
stars with log g < 3.5, defined as giants. Stars with Jc > 0.5 
are not plotted. Typical uncertainties are shown in the lower left- 
hand corners. 


dwarfs/sub-giants (log g ^ 3.5) and giants (log g < 3.5), 
which are shown in Figure 3b and 3c, respectively. 

5.1.1 Retrograde stars (J^/J^ < 0/ 

The stars with Jz/Jc < 0 (i.e. those moving retrograde with 
respect to the disc) are a-enhanced (([Mg/Fe]) ~ 0.4 — 0.5) 
up to a metallicity of ~ —1.5 dex. From there, their aver¬ 
age [Mg/Fe] ratios predominantly decrease with increasing 
metallicity, with stars reaching [Mg/Fe] 0.2 at [Fe/H] ~ 
—0.8. This is most clearly seen in Figures 3a and 3c. (There 
are also some stars with 0 < Jz/Jc < 0.5 that also follow 
this trend, albeit with a much larger spread.) 

The metallicity behaviour of our Jz/Jc < 0 stars is 
qualitatively similar to the trends found for surviving dwarf 
galaxies (e.g. Tolstoy et al. 2009). These stars also have 
abundances and kinematics that fit within our definition of 
accreted stars in R14 and §2, suggesting an accretion origin. 
Moreover, these trends are very similar to those found by 
Nissen & Schuster (2010), in which they investigate high- 
velocity stars in the solar neighbourhood. Finally, note that 
the majority of the stars in our sample with Jz/Jc < 0 sit at 
distances outside of the solar neighbourhood {Z > 1 kpc or 
R < 7 kpc or i? > 9 kpc). This is shown in Figure 4, where 
the positions of the stars are plotted in the Galactocentric 
R — Z plane. This further strengthens the conclusion that 
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Figure 4. Position of stars in the Galactocentric R — Z plane, highlighting stars with Jz/Jc < 0 as blue, triangles and those with 
0 J-l! 3c < 0.5 as red, circles, (a) All stars in our sample, (b) Only stars with log g ^ 3.5, defined as sub-giants/dwarfs, (c) Only 
stars with log g < 3.5, defined as giants. Typical uncertainties are less than 1 kpc in both R and Z. For the most distant stars, the 
uncertainties can range up to .^2 kpc. 


these are stars accreted from low mass dwarf gala:xies. It also 
proves that the results found by Nissen & Schuster (2010) 
are not exclusive to the solar neighbourhood. 


5.1.2 Low angular momentum stars (Jz/Jc < 0.5^ 

A second interesting trend can be identified at metallicities 
between —0.8 and —0.3 dex. Here, there appears to be a 
possible over-density of stars with Jz/Jc < 0.5 (both red 
and blue symbols in Figure 3) and [Mg/Fe] ~ 0.3 — 0.4. 
This is best seen in the giants (Figure 3c), showing also a 
slight decrease with increasing metallicity. This clump will 
be discussed further in §5.2. 


5.2 Identifying substructure in Jz/Jc 

In the previous section, we found stars that have both 
[Mg/Fe] ratios and specific angular momenta compatible 
with our chemo-dynamical analysis in §2. In addition, there 
may be an over-density of low- Jz/Jc stars at higher metal¬ 
licity and high-a. However, we have not yet determined the 
statistical significance of these features. Since they lie in the 
tail of the thin/thick disc of the Milky Way, they could be 
contaminated by thin/thick disc stars, or even completely 
spurious. 

To determine if these stars can be attributed to signifi¬ 
cant accreted components, and to account for uncertainties 
in their Jz/Jc values, in this section we investigate their 
detailed distributions in Jz/Jc^. This was initially done by 
slicing our data into a grid of 0.2 dex and 0.1 dex bins in 
[Fe/H] and [Mg/Fe], respectively. These were chosen so that 
the bin size in [Mg/Fe] was about twice the size of the un¬ 
certainty, while the bin size in [Fe/H] was chosen to retain 
significant number statistics (ristars > 10) in the majority 

® In addition to Jz/Jc, R14 also computed the vertical specific 
energy Ec) for each star. In that work, we found that no dis¬ 
tinguishing information could be extracted from the distribution 
in Bz/J'c (see, e.g.. Fig. 7 of R14). Thus, for the present analysis, 
we concentrate on Jz/Jc only. 
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Figure 5. Location of our stars in the [Mg/Fe] — [Fe/H] plane. 
The grid of bins used in the first step in the analysis in §5 is 
shown by the red, dashed lines. Note, at the extrema of the grid, 
the bins extend to infinity. 



of the metal-poor bins. Shifting the bins by 0.05 dex in any 
direction has no major impact on subsequent results. The 
grid is shown in Figure 5. At the extremes of both ranges, 
all stars are included to “infinity” in order to improve the 
statistics. 

For each [Mg/Fe] — [Fe/H] bin, we computed the com¬ 
plete distribution of Jz/Jc by snmming the PDFs of individ¬ 
ual stars present within that bin. The summed distribution 
was normalised such that the total area is equal to one. 
The grid of distributions for Jz/Jc for all abundance bins 
is shown in Figure 6. This figure reveals several interesting 
characteristics: 

(i) The majority of the [Mg/Fe] — [Fe/H] bins show a pri¬ 
mary peak at Jz/Jc ~ 1.0 with an extended tail towards low 
Jz/Jc values, i.e. the Galactic disc. 

(ii) The red histograms in panels Ic, 2c, 3c-d, and 4e indi¬ 
cate those bins with [Mg/Fe] < 0.3 (defined in §2) and that 
have a signature of a distinct low-Jz/Jc component in ad¬ 
dition to the peak at Jz/Jc ~ 1. The low- Jz/Jc component 
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Figure 6. The distribution of the vertical component of the specific orbital angular momentum (Jz/Jc) for each [Mg/Fe] — [Fe/H] bin 
defined in Figure 5. The distributions are computed by summing PDFs for Jz/Jc for all stars present within a given bin. This is then 
normalised such that the total area in a given bin equals one. The numbers in brackets give the total number of stars in each bin. The 
red distributions in Ic, 2c, 3c-d, and 4e show signatures of an accreted component that falls in line with our chemo-dynamical template 
defined in R14 (namely, [Mg/Fe] < 0.3 and Jz/Jc < 0.7). Additionally, panels 4a, 5a, 6b, 9d-f (shown in green) also contain possible 
low-Jz/Jc components. Panels which contain less than 10 stars, and thus suffer from low-number statistics, are shown in grey. 


peaks between 0.00 and 0.25 in all of the bins, except 4e. 
Panel 4e shows possible indications of a second component 
peaking between Jz/Jc ~ 0.2 and 0.5. 

(iii) There are also possible signatures of a low-Jz / Jc com¬ 
ponent in bins containing stars with high [Mg/Fe] ratios at 
metallicities above —0.8 dex (green panels). This is most 
evident in panels 4a, 5a, and 6b. Here, the stars are rela¬ 
tively metal-rich (—0.8 ^ [Fe/H] < —0.2) and a-enhanced 
([Mg/Fe] ^ 0.3). Panels 9d-f also contain a handful of low- 
Jz/Jc stars with high metallicities, [Fe/H] > 0.2. 

(iv) Bins that suffer from low-number statistics are shown 
as grey histograms. Although the statistics for these bins is 
poor, all of them contain low-Jz/Jc stars. 

(v) Panels la-b, 2a-b, and 3b all show low-Jz/Jc compo¬ 
nents. The [Mg/Fe] ratios of these bins are consistent with 
normal high-o? halo stars. 

The stars identified in panels Ic, 2c, and 3c-d contain 
the low-Jz/Jc stars at low metallicity hrst identified in §5.1. 
These stars are likely accreted stars, but their low metallicity 
and low Jz/Jc suggests they were likely accreted from a low- 
mass satellite galaxy (see discussion in §2). In contrast, the 


more prograde component in panel 4e is in line with the 
expected specific orbital angular momentum of accreted disc 
stars that have come from a higher-mass, low-inclination 
merger. We will discuss the statistical signihcance of this 
feature in §5.3. 

The panels 4a, 5a, and 6b confirm our detection of a 
low- Jz/Jc component at high-a and high metallicity in Fig¬ 
ure 3. These stars have [Mg/Fe] that are clearly distinct from 
those identified in (ii). The [Mg/Fe] ratios of these high-a 
stars show some similarity to stars observed in the inner disc 
(e.g. Bensby et al. 2011) and bulge (e.g. Bensby et al. 2013; 
Ness et al. 2013), and thus, may have similar origins. Fur¬ 
thermore, the stars in panels 9d-f could be stars that have 
migrated outwards from the inner parts of the Galaxy (e.g. 
Kordopatis et al. 2015). 

For the remainder of this paper, we will concentrate on 
identifying accreted stars and the accreted disc from point 
(ii). The low-Jz/Jc signal for high-a, metal-rich stars, de¬ 
scribed in point (iii), will be discussed in more detail in a 
forthcoming paper. 

Figure 6 shows the power of using the distribution of 
Jz/Jc in different abundance bins to pinpoint possible sub- 
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Figure 7. The same as Figure 6, but only stars that have probabilities, P(Jz/Jc ^ 0.7) < 0.05, are plotted. Panels 4d, 5d-e, 6d, and 
6f (now shown as magenta histograms) show a clear low-Jz/Jc component that was otherwise hidden in Figure 6. The number in the 
bracket gives the total number of stars in each bin. 


structure in a large sample of stars. These distributions pro¬ 
vide a much clearer picture than using only the abundance 
trends in Figure 3. 

Still, there could be an accreted component in the disc- 
dominated bins (blue bins with [Fe/H] (S —0.8 in Figure 6), 
which has been completely hidden by the low-angular mo¬ 
mentum tail of the disc component. In Figure 7, we cut 
away the Galactic disc stars by only plotting those stars 
which have a probability, P(Jz/Jc ^ 0.7) < 0.05. In other 
words, the stars plotted only have a five percent probability 
(according to their derived PDF) of having a Jz/Jc 0.7. 

Figure 7 reveals stars on low-Jz/Jc orbits in several 
disc-dominated bins at metallicities greater than —0.8 dex 
(magenta bins: 4d, 5d-e, 6d, and 6f). The peak value at 
Jz/Jc ~ 0 is relatively far away from our cut at Jz/Jc = 0.7. 
Thus, we can be confident that the peak is not just a result 
of this cut. These bins, however, contain very few stars and 
thus it is possible that some of these could be disc stars with 
very poor measurements. There are also several other bins 
(e.g. 7c-d and 8d) that also show indications of low-Jz/Jc 
stars. However, these bins are likely connected to the green 
histograms identified in Figure 6, which will be discussed in 
a later paper. 

It is clear from Figures 6 and 7 that some bins contain 


very few stars. This can make it difficult to assess the sig¬ 
nificance of any detected accreted component. In the next 
section, we combine bins to reinforce the statistics. 


5.3 Characterizing the accreted population 

In this section, we explore the evidence for accreted stars 
present in our sample. For this analysis, we only considered 
stars with [Fe/H] < —0.2. Above this metallicity, we do not 
expect to find a significant accreted component (see §5.2 as 
well as Fig. 1 of R14). In an effort to improve the statistics 
of the bins associated with accreted stars in Figure 6, we 
utilize a variant of our original chemical template in R14. 

In R14, we found that the abundance ratios of stars 
in surviving satellite galaxies reach an upper threshold of 
[Mg/Fe] ~ 0.3 in the metallicity range of the sample ana¬ 
lyzed in that work. We therefore used a cut at [Mg/Fe] = 0.3 
to define high- and low-a stars. However, a constant cut at 
[Mg/Fe] = 0.3 across all metallicities is not appropriate for 
our current Gaia-ESO sample that extends to higher metal¬ 
licity than R14. 

In Figure 8, we define three [Mg/Fe] — [Fe/H] bins: Bin 
A (-0.8 < [Fe/H] < -0.2), Bin B (-1.3 ^ [Fe/H] < -0.8), 
and Bin C ([Fe/H] < —1.3). The metallicity range of Bins 
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Figure 8. The sample plotted in the [Mg/Fe] — [Fe/H] plane, 
but now sliced into larger bins, similar to that defined in R14. 
Stars are separated into high- and low-a regimes, with a cut at 
[Mg/Fe] = 0.3 for stars with [Fe/H] < —0.8. At larger metallic- 
ities, the data were divided by a cut at [Mg/Fe] = 0.2, which 
was determined using Figure 7. The vertical, dashed lines di¬ 
vide our sample into three metallicity bins, shown as Bin A: 
—0.8 ^ [Fe/H] < —0.2, Bin B: —1.3 $5 [Fe/H] < —0.8, and Bin C: 
[Fe/H] < —1.3. These bins are also used in Figures 9—11. Stars 
with [Fe/H] > —0.2 (grey) are excluded from further analysis. 


B and C is similar to that in R14. We therefore applied the 
same cut at [Mg/Fe] = 0.3 used in R14 for these bins. For 
Bin A, this cut is clearly too high according to Figure 7 
where a low-Jz/Jc component is revealed for most bins with 
[Mg/Fe] < 0.2. We therefore adopted a cut at [Mg/Fe] = 0.2 
for this metallicity bin. 

We next computed the distribution of Jz/Jc for the 
high- and low-a stars (defined by the [Mg/Fe] cuts in each 
bin) in each corresponding [Mg/Fe] — [Fe/H] bin by sum¬ 
ming the Jz/Jc PDFs of the stars within each given bin. 
These are shown in Figure 9. Similarly to Figure 7, in Fig¬ 
ure 10 we explore the possibility that an accreted compo¬ 
nent could be hidden by the extended tail of the Galactic 
disc stars by plotting the Jz/Jc distributions for only stars 
with P(Jz/Jc ^ 0.7) < 0.05. This is most useful for Bin 
A, the high metallicity bin, which has the most substantial 
disc component. In the following subsections, we explain our 
findings. 


5.3.1 Bin A: -0.8 ^ [Fe/H] < -0.2 

In this bin, the high- and low-a distributions are very sim¬ 
ilar (Figure 9a) with a peak at Jz/Jc ~ 1 and a long tail 
extending down to Jz/Jc ~ 0. If we cut away the disc dis¬ 
tribution as shown in Figure 10a, we now see a peak around 
Jz/Jc ~ 0.1 — 0.2, which was hidden in Figure 9a. We also 
see a slight difference between the high-a and low-a compo¬ 
nents, with the low-a peak slightly shifted to lower Jz/Jc. 
A simple two-sample Kolmogorov-Smirnov test yields a p- 
value< 2 x 10“^®, implying that the difference is statistically 
significant. This suggests that even if the high-a component 
is just the tail of the disc population (spread to low Jz/Jc 
by uncertainties), the low-a component contains something 
else - a real low Jz / Jc component from an accreted satellite. 



Figure 9. Distributions of the vertical component of the specific 
orbital angular momentum (Jz/Jc) for the three [Fe/H] bins de¬ 
fined in Figure 8; (a) Bin A: —0.8 C [Fe/H] < —0.2, (b) Bin B: 
— 1.3 Si [Fe/H] < —0.8, (c) Bin C: [Fe/H] < —1.3. The distribu¬ 
tions are the sum of the PDF values within a given Jz/Jc bin, 
and are normalised such that the total area equals unity. The 
total number of stars whose PDF values were summed for each 
distribution is given in the legend. The high- and low-o distribu¬ 
tions, defined by the cut in [Mg/Fe] shown in Figure 8, are shown 
as black and red distributions, respectively. 


5.3.2 Bin B: -1.3 s/ [Fe/H[ < -0.8 

We find strong evidence for accreted stars in Bin B. Fig¬ 
ure 9b shows two distinct peaks in the low-a distribution: 
one with Jz/Jc ~ 0.8 and a low specific angular momentum 
peak at Jz/Jc ~ 0. This duality is not seen in the high-a 
distribution. Figure 10b further confirms what is seen in Fig¬ 
ure 9b, that there is a component in the low-a distribution, 
peaking at Jz/Jc ~ 0. 


5.3.3 Bin C: [Fe/H] < -1.3 

At the lowest metallicities, we see accreted stars that have 
low- Jz/Jc orbits, as shown in Figure 9c. Further, the Jz/Jc 
distributions shown in Figure 10c are almost unchanged, as 
neither has a disc-like component. At these metallicities, sur¬ 
viving dwarf galaxies show similar [Mg/Fe] ratios to stars in 
the MW (e.g. Tolstoy et al. 2009; Aden et al. 2011; Kirby 
et al. 2011). We might then expect that the accreted stars 
should also have similar abundances to the in situ popu¬ 
lation. We do see, however, subtle differences between the 
high- and low-a trends. For example, the low-a distribution 
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Figure 10. The same as Figure 9, but only stars with probabil¬ 
ities, P(J^,/Jc ^ 0.7) < 0.05, are plotted. The three [Fe/H] bins 
are defined in Figure 8; (a) Bin A: —0.8 ^ [Fe/H] < —0.2, (b) 
Bin B; -1.3 ^ [Fe/H] < -0.8, (c) Bin C: [Fe/H] < -1.3. The 
distributions are the sum of the PDF values within a given Jz/Jc 
bin, and are normalised such that the total area equals unity. The 
total number of stars whose PDF values were summed for each 
distribution is given in the legend. The high- and low-o distribu¬ 
tions, defined by the cut in [Mg/Fe] shown in Figure 8, are shown 
as black and red distributions, respectively. 

shows a peak at Jz/Jc ~ —0.5 in both Figure 9c and 10c, 
which is not seen in the high-a distribution. 

5 . 3.4 The age distributions 

In Figure 11, we plot the age distributions of both the high- 
and low-a stars with P(Jz/Jc ^ 0.7) < 0.05. As with the 
Jz/Jc distributions in Figure 10, these distributions are com¬ 
puted by summing the age PDFs of the stars within a given 
[Mg/Fe] — [Fe/H] bin. There are far fewer stars used to com¬ 
pute the age distributions, since we were limited to only 
those stars with log g > 3.5 (see §4.2). Even though we are 
strongly affected by low-number statistics, the distributions 
are still useful as consistency checks of our Hndings in the 
Jz/Jc distributions. 

The low-a stars in the low metallicity bins - Bin B 
(Figure 11b) and Bin C (Figure 11c) - show predominantly 
old ages (> 10 Gyr), consistent with a population of ancient 
accreted satellites. In contrast, the low-a distribution in the 
more metal rich Bin A has a peak age at ~ 6 — 9 Gyr, albeit 
with quite a large spread. Since lower-mass satellite galaxies 
are expected to grow in metallicity more slowly than the 
MW, we might expect that metal-rich accreted stars should 


Figure 11. Distribution of the ages of stars with log g IS 3.5 and 
probabilities, P{J^.IJc ^ 0.7) < 0.05, in the three [Fe/H] bins 
defined in Figure 8; (a) Bin A: —0.8 ^ [Fe/H] < —0.2, (b) Bin B: 
— 1.3 Si [Fe/H] < —0.8, (c) Bin C: [Fe/H] < —1.3. The total num¬ 
ber of stars whose PDF values were summed for each distribution 
is given in the legend. The high- and low-o distributions, defined 
by the cut in [Mg/Fe] shown in Figure 8, are shown as black and 
red distributions, respectively. 

have younger ages than their more metal-poor counterparts. 
Thus, it is possible that some of these stars were accreted. 
However, such ages are also typical of normal thin and thick 
disc stars (e.g. Haywood et al. 2013; Bensby et al. 2014). 
These can extend to low alpha - especially if the stars are 
members of the thin disc or have migrated inwards from 
the outer disc (e.g. Sellwood & Binney 2002). It is therefore 
probable that some of the low-a stars in Figure 11a are 
actual disc stars. This means that even with our probability 
cuts, this bin still suffers from disc contamination. 

5.3.5 The accreted component 

At metallicities less than —0.8 dex (Bins B and C), we have 
identihed a clear signature of an accreted component. The 
metallicity range and peak Jz/Jc ~ 0 of this component, 
suggests that the stars were either accreted from one or 
more low-mass satellite galaxies or one or more LMG-like 
mass satellites on very high-inclination orbits (see §2, as well 
as R14). This latter interpretation seems unlikely, however, 
since such a merger ought to deposit also higher metallicity 
stars ([Fe/H] > —0.8) at low [Mg/Fe] with low Jz/Jc. The 
Jz/Jc distribution in Figure 10a (Bin A) shows a higher 
density for Jz/Jc > 0.2 than that shown in Bins B and 
C. Further, most of the stars in Bin A have ages similar 
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to Galactic disc stars. It is therefore likely that the distribu¬ 
tions in Bin A are mostly comprised of disc stars and contain 
few accreted stars. We conclude that although there may be 
some accreted stars in Bin A, we do not find a significant ac¬ 
creted component in this metallicity bin, implying that the 
accreted stars in Bins B and C came from low-mass satellite 
gala:xies. 


6 MEAN MASS OF THE MERGING 

SATELLITE(S) 

In §5, we confirmed that our sample contains accreted stars. 
We suggested that these stars were likely accreted into the 
halo of the MW from low-mass satellite galaxies. In this 
section, we attempt to estimate the mass of the merging 
satellite galaxies from which the accreted component in our 
Gaia-ESO sample originated. 

In Kirby et al. (2013), the authors compiled the metal¬ 
licity of individual stars in several dwarf irregular and dwarf 
spheroidal galaxies in the Local Group. From this data, they 
produced a stellar mass—metallicity relation for the Local 
Group dwarf galaxies. In addition, they combined this with 
the mass-metallicity relation for more massive galaxies ob¬ 
served in the Sloan Digital Sky Survey (Gallazzi et al. 2005). 
We adopted this relation to estimate the mass of the merg¬ 
ing satellites. We note, however, that the Kirby et al. (2013) 
mass—metallicity relation only depends on the mean metal¬ 
licity of a dwarf galaxy. They do not consider the full range 
of metallicities of the dwarf galaxy, which could provide ad¬ 
ditional Information about the mass distribution and star 
formation history of the satellite. 

In §5.3, we showed that there is an accreted, low- 
a ([Mg/Fe] < 0.3) component at metallicities [Fe/H] < 
—0.8. In addition, there could be some accreted stars with 
[Mg/Fe] < 0.2 at metallicities, —0.8 ^ [Fe/H] < —0.2. If 
we now consider all of the low-a stars at these metallicities 
that have P(Jz/Jc ^ 0.7) < 0.05, the average metallicity 
is ([Fe/H]) ~ —1.0. Using the mass-metallicity relation in 
Kirby et al. (2013), this average metallicity corresponds to 
a satellite galaxy with a stellar mass of M* ~ 10®’^ Mq. If, 
on the other hand, we assume that the stars with metal¬ 
licities [Fe/H] > —0.8 are all disc stars, then the average 
metallicity decreases to ([Fe/H]) ~ —1.3. This corresponds 
to a satellite galaxy with a stellar mass of M* ~ lO^ "* M 0 . 

Both of these values for the stellar mass are substan¬ 
tially less than the LMC, which has a stellar mass of 
M,,lmc 10®'® MQ(van der Marel et al. 2002). These satel¬ 
lite galaxies would therefore suffer very little disc plane drag¬ 
ging and we would expect them - when averaging over sev¬ 
eral mergers - to deposit their stars on nearly non-rotating 
orbits. It is reassuring, then, that the Jz/Jc distribution for 
these accreted stars in our sample peaks at zero. 

It is difficult to quantify exactly how many mergers our 
accreted star signal corresponds to for two reasons. Firstly, 
there is a degeneracy between the number and stellar mass 
of the mergers that is only weakly broken by the metallic¬ 
ity information (see discussion above). Secondly, to properly 
use information from the Gaia-ESO survey about the num¬ 
ber density of stars, we must correct for the survey selection 
function. This requires a good understanding of the stellar 
populations both in the Milky Way and the accreted satel¬ 


lites. Such detailed modelling is beyond the scope of this 
present work. 


7 THE QUIESCENT MILKY WAY 

In §5 and §6, we confirmed our previous results from R14. 
A key difference from our previous work, however, is that 
we do now clearly detect low-mass accretion events. That 
our analysis finds these, but not any at higher mass proves 
that the MW has suffered few if any LMC-mass mergers 
since its disc formed ~ 9 Gyr ago. As a result, it will not 
have a significant dark matter disc. This is reinforced by 
our simulation results in Figure 1 where we see that using 
just the kinematics alone, a massive LMC-like merger would 
stand out as a distinct peak in Jz/Jc space that we do not 
see in the Gaia-ESO data. 

Our result that the MW has had no significant mergers 
since its disc formed is supported by many complementary, 
albeit more indirect or model dependent, studies of the MW. 
Hammer et al. (2007) investigated the stellar mass and an¬ 
gular momentum of the MW, as well as the [Fe/H] of stars in 
its outer regions, finding that it did not undergo any signif¬ 
icant mergers over the last ~ 10 Gyr. The kinematic groups 
identified by Helmi et al. (2006) that are likely formed from 
the accretion of stars from satellites also support this notion, 
since they have been shown to consist of predominantly old 
stars (e.g. Stonkute et al. 2013; Zenoviene et al. 2014). Re¬ 
cently, Casagrande et al. (2015, submitted) used the com¬ 
bination of asteroseismology and spectroscopy to study the 
vertical age gradient in the Galactic disc. They found no 
signatures of over-densities in the age distribution for ages 
younger than 10 Gyr, which is consistent with a quiescent 
MW. In addition, Minchev et al. (2014) found that the veloc¬ 
ity dispersions of stars in the Galactic disc decrease with in¬ 
creasing [Mg/Fe], suggesting that the MW has experienced 
few late major mergers. This has recently been confirmed 
with the Gaia-ESO Survey (Guiglion et al., submitted). Fi¬ 
nally, the nearly spherical shape of the dark matter halo, 
inferred from measurements of the local dark matter den¬ 
sity (e.g. Read 2014, and references therein), and the steep 
fall-off in the stellar halo density profile (e.g. Deason et al. 

2013) both suggest that the MW had a quiet merger history 
with an insignificant dark matter disc. 

In contrast to the above studies, two over-densities - 
the Monoceros Ring and Triangulus Andromeda - have 
been identified in the disc. These contain intermediate-age 
stars (e.g. Carraro 2015, and references therein) and it has 
been suggested that they could result from accreted satel¬ 
lite galaxies (e.g. Gonn et al. 2012). If true, this would imply 
more recent mergers, at odds with our results here. However, 
both features could also be the result of disc oscillations (e.g. 
Xu et al. 2015) or flaring (e.g. Lopez-Corredoira & Molgo 

2014) . Our results here disfavour a merger origin for these 
features and lend support to the disc oscillation/flaring sce¬ 
nario. 


8 GONCLUSIONS AND FUTURE PROSPEGTS 

In this work, we expanded the chemo-dynamical analy¬ 
sis developed in Ruchti et al. (2014) to a large sample of 
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4,675 MW stars in the third internal data release from 
the Gaia-ESO Survey. Our analysis uncovered a significant 
component of accreted stars on nearly non-rotating orbits, 
which were likely accreted from satellite galaxies with stellar 
masses of — 10®'^ M©. However, we found no evidence 
of a significant prograde accreted disc component, which 
would be the result of one or more ^ LMC-mass mergers. 
We conclude that the MW has no significant accreted stel¬ 
lar disc, and thus possesses no significant dark matter disc. 
This implies that the MW did not experience any substantial 
( 1 : 10) mergers after the formation of the Galactic disc 

9 Gyr ago), and therefore has had a relatively quiescent 
merger history. 

On the other hand, our template revealed a possible 
signature of a component of metal-rich (—0.8 ^ [Fe/H] < 
—0.2), Q-enhanced ([Mg/Fe] ^ 0.3) stars with low specific 
vertical angular momentum. In a follow-up paper, we will 
investigate the nature of these interesting stars using com¬ 
parisons to recent state-of-the-art cosmological simulations 
(e.g. Agertz & Kravtsov 2014). 

Finally, our results here have interesting implications for 
the local dark matter density pdm and the local dark mat¬ 
ter velocity distribution function /dm- Firstly, the lack of a 
dark disc suggests that /dm is near-Gaussian (rather than 
the double-Gaussian distribution expected in the case of a 
significant dark disc component; e.g. Read 2014). Secondly, 
if Pdm is found to be substantially higher than spherical ex¬ 
trapolations from the rotation curve would imply, the lack 
of a rotating dark disc allows us to say with confidence that 
the dark halo of the MW is locally oblate (e.g. Read 2014). 
Such a locally oblate halo is expected due to contraction of 
the dark halo in response to disc growth (e.g. Katz & Gunn 
1991; Debattista et al. 2008; PifH et al. 2015), but evidence 
for it remains elusive (e.g. Read 2014). 
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